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We report theoretical and experimental results on 4-wave capillary wave turbulence. A system 
consisting of two inmiscible and incompressible fluids of the same density can be written in a 
Hamiltonian way for the conjugated pair (?7,^'). When given the symmetry z — > —z, the set of 
weakly non-linear interacting waves display a Kolmogorov-Zakharov (KZ) spectrum ~ fc"* in 
wave vector space. The wave system was studied experimentally with two inmiscible fluids of 
almost equal densities (water and silicon oil) where the capillary surface waves are excited by a 
low frequency random forcing. The power spectral density (PSD) and probability density function 
(PDF) of the local wave amplitude are studied. Both theoretical and experimental results are in 
fairly good agreement with each other. 
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^SJ , Introduction: Wave turbulence []| deals with a set of nonlinear random waves in a dispersive medium that, although 
forced far from thermodynamic equilibrium, can be described statistically. This description is done by means of a 
kinetic equation for the spectral density distribution Uk which evolves through resonant interactions of Af waves. 
Besides equilibrium solutions represented by Rayleigh-Jeans distributions, the kinetic equation can display stationary 
^ [ power-law nonequilibrium solutions rik ~ k^^, i-i > 0, called Kolmogorov-Zakharov (KZ) spectra describing the energy 
exchange (or other conserved quantities) between large and small scales. KZ spectra has been predicted theoretically 
and observed numerically and experimentally in systems such as bending waves in elastic sheets 2, 3], Alfvcn waves 
in plasmaQ, spin waves in solidsQ and gravity waves in fluids 0, 0, H 0|, to name a few examples. In all the above 
, systems, although the theoretical description depends on several strong constraints (negligable viscosities, density 
(~| ■ constrasts, aspect ratios, etcP, [l3|), experimental and numerical results corroborate the theoretical prediction of 
O . the appearence of power-law nonequilibrium spectra. Still, there are certain features that are yet to be studied and 
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compared between theory, experiments and numerics, for instance, the non-gaussianity of the wave amplitudes 14 1. 
the nature and existence of intermitency in a wave svstem[l5l.[l6j. the role of symmetries and dissipation in the wave 
interactions 17, l3| or the deviations of the exponent u from the theoretical value. 
I In this Letter we focus on capillary wave turbulence 13, 11, l3l and the effect of symmetries in the wave interactions. 
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We study theoretically and experimentally the statistical properties of random waves at the interface between two 
inmiscible and incompresible deep fluids of equal densities (pi and P2) and depths {hi and /i2). Due to these facts 
^f) • the symmetry z — > — ^ is forced on the system: the typical 3- wave capillary wave turbulence breaks down and a 
four-wave resonant interaction appears as the leading order perturbation. We discuss the effect of this symmetry on 
' the nonlinear type of wave interaction and on the KZ spectrum of the wave amplitude. We compare the theoretical 
^\ ] prediction with the experimental measurement of the power spectral density (PSD) of the local wave amplitude at 
' a water-oil interface in the limit of the Atwood number A — (pi — p2)/{pi + P2) 0. Also, the probability density 
function (PDF) for the excited surface wave amplitude is computed. Both results are contrasted with gravity-capillary 
, , ■ wave turbulence measurements. The level of agreement between theoretical and experimental results stresses the fact 
rS ] that capillary wave turbulence is a robust phenomenon for nonlinear random waves. 

' Theoretical study: Let us study the system of potential flows of two incompressible and immiscible fluids in a box 
of height 2h = hi + h2, where r](r, t), r = [x, y) corresponds to the surface elevation between them, pi is the density 
of the bottom fluid {—hi < z < rf), p2 the density of the upper fluid (77 < z < ft.2) with pi > p2 and a the surface 
tension coeflacient between the two fluids. The flows are defined by the velocity potentials 0i(r, z, t) in the lower fluid 
and (f)2{j',z,t) in the upper fluid with Vi^i = wi,V02 = V2- It is possible to prove that the dynamics of the two fluid 
interface have a Hamiltonian structure [Til [l9| . i.e., 

dr^ir) _ 5H 9«'(r) _ 5H 
dt ^ dt ^' 
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with ^'(r) = pi(f)i{r, r]{r)) — P202(?', '7('')) a-nd H ^ K + U given by 
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Pi ^ dzdr 

-hi ^ 



fl2 dzdr 
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dr 
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where 1. correspond to the r coordinates. It is easy to see that the Hamihonian corresponds precisely to the sum 
of the kinetic energy K (upper line) and potential energy IJ (lower line) of the system. The system is also con- 
strained to the bondary conditions dzi>\\z=-hx = dz4>2\z=h2 — (zero normal velocity at the container walls), 
[(■y2_L V_L)r7(r) — V2z\z=ri ~ [(■f^i-L V_L)?7(r) — viz]^^^ (continuity of the normal velocity at the interface) and the in- 
compressibility conditions V^^i = V^(/)2 = 0. A formal expression can be found for the Hamiltonian [S^. The 
original work was presented for the 2D case, but can be easily extended for the 3D case. The kinetic energy can be 
expressed as 



K = / *G2(??)(p2Gi(7?) + piG2(^))-'Gi(77)^dr, 
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where Gi (77), 62(77) are the Dirichelet-Ncumann operators for the fluid domain ^hi < z < ri[r) and ry(r) < z < h2 

respectively, defined by Gi{r])(j>i{r,rj{r)) = (—1)' Vr(/'iVr??(r) — ^ . In 3D systems it docs not seem possible to 

write an explicit Hamiltonian in terms of rj and ^P. This problem is bypassed by using the "small angle aproximation" 
to write the Hamiltonian as an infinite Fourier series in fc-space[ll [3, [llj. In terms of the operators Gi(r]),G2ir]), 
it corresponds to find an asymptotic series in term of the small parameter krjQ -C 1, with k the wave vector and 
r]Q the characteristic surface elevation. When adding the symmetry z —z to the initial problem (in this case by 
impossing equal depth and density) the expansion naturally needs to satisfy this constrain. Therefore, the order of 
the nonlinearity increases from A/'=3 to 4, and also the system becomes gravity free. In the Hamiltonian expansion 
H = H2 + Hi..., one gets 



H2 = 
Hi = 



1 



kta.VLh[kh]'l! k'^ _k + crk^VkV- 



dk 



J (^l,2;3,4*fei ^'fe2 Vk3 Vki + T^%3,4:Vki ?7fc2 Vks Vk, 

X 5{ki + k2 + ks + ki)dki2zi, 
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where = 0, which comes explicitly from the the Dirichelet-Neumann operators expansion [20j. Physically, the 
z —z symmetry is not merely imposed by eliminating gravity in the limit A — > 0, but also by impossing equal 
depth of both fluids, making the capillary surface waves unable to disinguish up from down in the vertical direction. 
Even more, using the same symmetry arguments, every odd term in the expansion of H will be zero in this limit. 
Therefore, the computed transfer matrixes 

rS;3,4 = Y^(-2tanh[fci/i] ((fci • fc2)/ci) - 2tanh[fc2/i] ((fe • ^1)^2) 



\_ / [fci.(fc4 + fcl)][fc2-(fc4 + fcl)] 
16p V |fcl + fc4 

\_ / [fci-(fc3 + fcl)][fc2-(fc3 + fci)] 

16p V Ifci + fcal 

^1,2:3,4 = --^(fcl ■ fe)(^3 ■ ^4), 



cotli[|fci + ki\h\ 
cotli[|fci + fcaj/i] 



[fc2-(fc4 + fc2)][fcl-(fc4 + fc2) 
|fc2 +fc4| 

[fe ■ (fcs + A:2)][fci • (fca + /C2) 



■ coth[|/s2 + ki\h\ 
■coth[|A;2 + kn\h\ 



(5) 



depend on four wave vectors. 

Following [l|, 01, we write ??A:=^ X^s ^fc ^^"^ ^''^ ^^lo in canonical variables A\. where s=± such that 

= Ak, A~^ = A*_f. and Xk~{ ^'^2^]^k^ ^ . From the hamiltonian evolution of A|, a hierarchy of linear equations 
for the averaged moments ( (A^M^^^), (A^^^ A^^) and so forth) is written. An assymptotic closure can be given 
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when there exists a separation of Hnear and non-hnear time scales and the system can be regarded as homogeneous 
in space[l^. The wave spectrum Uk, that comes from the second order moment (A^j A^^) = Uk^S'-^^^ki + k^) satisfies 
a Boltzmann-type kinetic equation describing the slow evolution of the wave spectrum through a four-wave resonant 
process: 



d 



Sl,S2,S3 

1111 

h Si h 32 h S3 

nk rik^ 
X S{ujk + siujki + S2UJk2 + satoka) 

X 6ik + ki+ k2 + k3)dki23 (6) 

where iIfe'/l';_l3ffc^=-isi^Pi234[siS2^^T^^ 2^^^ is the scattering matrix, 7'i234[-] is the 

sum over the twelve possible permutations of 1, 2, 3 & 4 and e is a small parameter related to weak non-linearities which 
makes possible the time-scale separation. With this kinetic equation, we seek isotropic non-equilibrium distribution 



solutions 2l| . Despite some differences with the usual kinetic equation, the method of Zakharov can be applied here as 



0. In the deep fluid limit, one finds that the scattering matrix L and frecuency tOk ~ ^ / j-fc"^ tanh[fcft,] of capillary 



2p' 

waves are homogenous functions of degree (3 = 3 and a = 3/2 respectively, i.e. ^Afc Afei'Afc2 Afea ^ '^''^fc fcV fcl 'fca ■ 
Looking for a power-law solution of the form nk = Ak~^ with A a constant, it is possible to perform the Zakharov 



transformation over the right hand side of Eq.®, called coUisional terms [2lJ. In such way one gets a stationary out 
of equilibrium spectrum with /i = 4, that represent a constant energy flux solution. If wc consider a flux of energy 
per unit of mass P (it has dimensions of (length/time)^ through the big scales towards the small scales, one can 
find an explicit expression for A that leads to Uk = CP^/^ pk~^ . C is a pure real number that depends on some 
integrals directly related with the coUisional term that can, in principle, be computed numerically. No inverse cascade 
is allowed in this system due to the 3 —s- 1 wave interaction process, as it was already reported in Q. As wc are 
considering a statistically homogenous system in space, it is natural to compute the moments by taking space averages. 
Nevertheless, from the experimental point of view, taking space averages is quite a difficult task. On the contrary, 
time averages of local properties are much more accessible. It is possible to relate both for stationary solutions in 
the linear regime as nT{uj)duj oc n{k)k'^~^dk, where nT{oj) is the time averaged wave number, in frecuency domain. 
Using LUi; of capillary waves one gets nxiu)) oc P^^^uj~'^^^ . Thus we obtain the spectrum for local surface elevation 

(hc.|')TOcPl/3^-8/3. 

Experimental study: A Plexiglass container (height h=60 mm, lenght ^=100 mm, depth d=80 mm) is half filled 
with destiled water (density pi =1.00 g/cm^, kinematic viscosity j^i =0.01cm-^/s) and half filled with silicon oil 
(PDMS AB112153 from ABCR, density p2 =0.93 g/cm^^, kinematic viscosity 1^2 =0.07cm^/s). The intcrfacial tension 
coefficient cr ~ 30 mN/m[22l . The equilibrium interface position is measured at 35 mm. The capillary surface waves arc 
excited by a wave-maker that plunged completely into the upper fluid, oscillating vertically. The wave-maker is driven 
by an electromagnetic vibration exciter via a power amplifier. The random forcing, supplied by the source output 
of a dynamical spectrum analyser, is low-pass filtered between 0-fdriv=3 Hz. The excited surface wave amplitude 77 
is locally measured 4 cm away from the container walls by means of a wire capacitive gauge, 0.1 mm in diameter. 
The measured capacitive fluctuations are proportional to the local wave amplitude ones. They are sampled at 800 
Hz during 300 s, and filtered numerically at 500 Hz to avoid aliasing. The only noteworthy difference in this set-up 
that in this case both the dielectrics are liquids of similar densities and similar viscosities. We have checked the linear 
response in the local wave amplitude rj by changing fluid depths and the constant frequency response (in magnitude) 
of the wire probe for the water-oil boundary in a frequency band between 1 to 100 Hz. 

With the acquired data, the probability distribution function (PDF) of ry, normalized by its rms fluctuations tr,, 
is calculated, as shown in Fig. [T] (main). Notice that (77) and that its fluctuations arc close to being symmetric 
with respect to 7/ = 0. No exponential tails arc found. The kurtosis is slightly larger than 3, but not large enough 
to exclude gaussianity. For comparison, wc show in Fig. [T] (inset) the PDF of r^/cr,, when gravity-capillary wave 
turbulence developes. We see a clear asymmetric tail (positive skweness) as it has been shown elsewhere This 
contrast is a clear indication of the symmetry imposed in the system: there is no external field (such as gravity) 
that breaks the z — > — z parity so the surface perturbations are symmetric with respect to 77 = 0. It is unclear if 
the fluctuations are indeed gaussian, but resolution of large events could not be made in the present experimental 
set-up. The wave system has a very low Atwood number A = (pi — P2)/ (pi + P2) —0.04, reducing the effective gravity 
drastically. One finds that in this system the capillary length Ic = 2Try^ Aa/ g{pi + P2) where the crossover from 
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FIG. 1: (Color online) Probability density function (PDF) of the normalized local wave amplitude ry/cr,, at the interface of two 
inmiscible fluids with A ~0.04 (full blue line) and a parabolic fit (dashed red line). Inset; PDF of rj/a^ at the interface of a 
water-air interface for A=l (full blue line) and a parabolic fit (dashed red line). 

gravity to capillary regime takes place is an order of magnitude larger than in a liquid-air interface problem psj. The 
frequency crossover between gravity and capillary regimes fc ~ lOc/It^ ~ tt-^/ Ag/2lc will be obtained at a frequency 
close to 3-4 Hz. Therefore, when the frequency cut-off of the forcing is larger than fc, the only KZ-type spectrum we 
can observe is the capillary one. In Fig. [2] we show both the pure capillary (main figure) and the gravity-capillary 
(inset) powerlike spectra. For pure capillary waves, as the forcing amplitude is increased, low frequency normal modes 
and harmonics of fdriv dissapear and a power-law spectrum developes. Only one scale-invariant spectrum with slope 
^ —2.75 ± 0.05 appears in the capillary-driven transparency window (for frequencies larger than the characteristic 
frequencies of the broad-band forcing), which is within the experimental error from the theoretical f~^^^. In this 
wave turbulent regime, no cusps over the wave crests were observed, which sustains the assumption krjQ and 
eliminates the possibility of singularities polluting the spectral content of the signal. The gravity-capillary spectra of 
Fig. [^l^inset) is calculated from the local elevation of surface waves when the lighter fluid (PDMS) is removed. Waves 
are excited using a wavemaker pluging in water driven by a low- frequency {fdriv=i Hz) random forcing similar to 
[^, [l^. Both gravity /-5.35 capillary /^^'^^ wave spectra are within the experimental range of 0]. 

Conclusions: In this paper we show the appearence of 4- wave capillary turbulence at the interface of two inmisci- 
ble and incompresible fluids with almost equal densities, due to an imposed spatial simmetry. The theoretical wave 
spectrum of the amplitude fluctuations behaves as f~^^'^ in frequency domain and fc~'* in wave vector space. Experi- 
mentally, we have computed the PDF and PSD of the local amplitude fluctuations 77 at the interface of an oil-water 
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Frequency [Hz] 

FIG. 2: (Color online) Power spectral densities (PSD) of the local wave amplitude rj at the interface of two inmiscible fluids 
with A ~0.04 (full blue line) for low (bottom) and high (top) forcing amplitudes. Best fit slope -2.75 (dashed red line). Inset; 
PSD of 77 at a water-air interface for A=l (full blue line) and best fit KZ spectra (dashed red line) for gravity (-5.35) and 
capillary (-2.52) waves. 

mixture. The computed PDF of 77 has a gaussian form, and no exponential tails where found0]. The experimental 
PSD of 77 shows a power-law behavior ^ /'^ ''^ in agreement with the expected theoretical slope for 4- wave interaction 
process. The excited surface waves in the experimental set-up display wave turbulence-type behavior which agrees 
with the theoretical prediction, even though A ^0 and both viscosities are not small, as the theory necesitates. A 
point worth studying both experimentally and theoretically is their effect on the dissipation scale and the effect of 
different depths, specially in the limit /ii//i2 1. It must be noticed that the computed spectrum is close to the 
theoretical /""/^ for 3-wave interaction processes. This slight departure can be explained by the fact that the 3-wave 
process may not suppressed completely due to a slight asymmetry in densities. The study of the slope dependence on 
A will be presented elsewhere. 
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